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Abstract. A study of the use of fractals in top non-leptonic decays for the sake of discrimination against background is 
presented. Preliminary results show that fractals may provide a useful check for top event enrichment techniques. 
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1 Introduction 

The lack of special features in top quark non-leptonic decays makes the extraction of this signal particularly hard to 
tackle with conventional techniques. The exploitation of decay particle correlations provides a possible way out to solve 
the problem. This is the case of neural-network classifiers, which have been used with some success [1-3]. Another tool 
may be provided by fractal geometry, which has already been applied to particle physics in other contexts, like e.g. in 
[4]. 

In this preliminary study we explore the capabilities of fractal geometry concepts in helping to discriminate 
between top quark non-leptonic decays and the underlying background. For this sake we concentrate on the 1 1 -> non- 
leptonic channel inpp interactions at the energy of the Fermilab Tevatron. We look for fractal quantities which exhibit 
the most pronounced differences when calculated on 1 1 -> non-leptonic and on background events. It turns out that for 
some of these quantities such differences are appreciable. They are not large enough to establish and identify by 
themselves a top signal, however they may provide useful checks for results obtained by other techniques. 

2 Fractal dimensions 

Fractal geometry is widely used in a variety of problems and does not need an introduction (for a recent particle physics 
application see, e.g., [4]). We only define entities which will be used in the present study. They are related to the 
concept of fractal dimension. In order to define it let us first consider simple examples based on uniform mass 
distributions. 

If we consider a mass distribution over a line starting at R = 0, the mass M(R) of the line segment up to R has a 
behavior 

M(R)ocR' (1) 

For a mass distribution uniform on a plane, one similarly has 

M(R) oc R 2 (2) 

where R now is the radius of a circle. 

The exponent of R takes the name of fractal dimension. If we take the massive line and we consider it on a plane, 
the dependency M(R) oc R 1 still holds with R now representing the radius of a circle centered on the origin of the line. 
Thus its fractal dimension provides a geometrical characterization of the system. 

For a mass distribution intermediary between the two considered, if it turns out that 

M(R) oc R D (3) 

D having some non-integer value, we are still in the presence of a scaling law characteristic of a pure fractal structure. If 
we plot log M(R) versus log R we obtain a straight line whose slope is D. 

If D=D(R), i.e. D depends on R, the distribution is said to represent a multifractal. Although a scaling law does 
not hold anymore, it may be still useful in such a case to employ the entities which are introduced for the handling of 
fractals. 

In a particle physics application the mass is to be substituted with something else associated with the particles. In 
this study we consider particle weights given by the particle transverse energy E Ti or a power of it E T f . As to R we 
define it as the particle pseudorapidity r\ or as the particle azimuth <j). More exactly we take 

R = Climax (4) 

or 

R = Wk (5) 
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where r\ max is the maximum possible r\, according to experimental acceptance, and -% < <p < Ttis measured with respect 

— > 
to the minor event axis m defined in the following section. 

We define the quantity 

y E q 

I{R^= ^< R a Tt (6) 

I— I all ' * 

R in the numerator represents a lower cutoff which may turn out to be useful so as to exclude particle contributions 
from some r\ or <p intervals. If R > it is in general I < 1. 

One may define a R dependent fractal dimension by 

q dlogR 
The differentiation, however, leads to higher statistical errors. In the applications we consider in this study it is thus 
better to directly use the quantities I(R,q) of eq. 6. 

3 Event selection 

tt and QCD background events are provided by the event generator COJETS 6.24 [5]. For pp interactions at a cm. 
energy E cm =1.8 TeV and a top quark mass m t = 175 GeV it yields a total tt cross-section of tt ) = 7.31 pb. 

A series of cuts is imposed in order to: i) suppress the contribution of leptonic and semi-leptonic decays of the 
top quark; ii) reduce the amount of the background. Where possible, the cuts mimic those used in the neural-network DO 
analysis [1-3]. 

For the suppression of the top leptonic and semi-leptonic decay contributions, events containing an electron or 
muon with a transverse momentum p T > 20 GeV or a missing transverse energy Ej lss > 20 GeV are excluded. 

For the sake of background reduction the following requirements are imposed: 

i) Total transverse energy Ej > 120 GeV; 

ii) At least one muon with transverse momentum p T > 4 GeV and absolute pseudorapidity I r\ \ < 2.5. 

iii) A linear aplanarity A L > 0.5. Considering only particles with 1 77I < 2.5 (neutrinos excluded), after finding the 

— > — » 

thrust, or major axis M, in the transverse energy plane, and defining the minor axis m as the axis orthogonal to it, linear 

— > — > — > — » — » 

aplanarity is defined as A L = £j | E t -m \ I Y,i I E i ■ M |, where E , is the transverse energy vector of particle i . 

For the integrated luminosity L = 1 10.2 pb" 1 applying to the DO event sample, the events left after all the cuts are: 

i) 128 events for tt ; ii) 3868 events for the background. 

For the present study we employ 1000 events for tt and 1000 events for the background. 

4 Fractal study for top non-leptonic decays 

We have looked for I(R,q)'s which exhibit the most pronounced differences between tt -> non-leptonic and background 
events. We have considered R 's defined by eq.s 4 and 5. In eq. 6 we have set log R = -Vi, which eliminates r\ and <p 
regions where I(R,q)'s are negligible. A range of integer values has been considered for q. 

5 Results 

Fig.s 1-4 show a selection of the results of this preliminary study. From them it appears that some of the fractal 
quantities I(R,q)'s exhibit appreciable differences between tt -> non-leptonic and background events. 

Comparisons between fig.s 1,3 and 2,4 show that differences are more pronounced for rather than for r\ . 



R. Odorico: Fractals in top non-leptonic decays 

Comparisons between fig.s 1,2 and 3,4 show that differences become considerably higher for high £V particles. 

5 Conclusions 

From this preliminary study, one can conclude that fractal differences between 1 1 -> non-leptonic and background 
events may help to discriminate between the two channels. They do not appear to be strong enough to allow by 
themselves to establish or to identify a signal, however they may represent a useful check for other signal enrichment 
techniques like, e.g., neural networks [1-3]. 
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Fig 1. Comparison of I(R,q=l)'s, with R defined as in eq. 4, for tt — > non-leptonic and 
background channels. 
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Fig 2. Comparison of I(R,q=l)'s, with R defined as in eq. 5, for tt — > non-leptonic and 
background channels. 
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Fig 3. Comparison of I(R,q=3)'s, with R defined as in eq. 4, for tt — > non-leptonic and 
background channels. 
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Fig 4. Comparison of I(R,q=3)'s, with R defined as in eq. 5, for tt — > non-leptonic and 
background channels. 



